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Abstract 

We present a mean field analysis of a fully frustrated Ising spin 
model on an Ising lattice gas. This is equivalent to a degenerate 
Blume-Emery-Griffiths model with frustration, which we analyze for 
different values of the quadrupolar interaction. This model might be 
useful in the study of structural glasses and related systems with 
disorder. 

PACS: 05.50 



1 Introduction 

In the last two decades the physics of complex systems, ranging from dilute 
magnets to structural glasses has been captured by models which couple Ising 
variables with lattice gas or Potts variables [fl-IMI) i-e. models with this type 
of Hamiltonian: 

- j3H = ^ JeijSiSjniUj + ^ Krurij + /j ^ nj, (1) 
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where Eij = ±1 are quenched variables associated to pairs of nearest 
neighbour sites, J > is the interaction between the Ising spin variables {S = 
±1), K is the interaction between the particles, = 0, 1 are the lattice gas 
variables, n is the chemical potential. The spins can interact each other 
ferromagnetically (e^ = 1) or antiferromagnetically {Sij = —1). 

For Eij = 1 everywhere, this model concides with the original Blume- 
Emery-Griffiths model (BEG) 0-[|l3] with an extra degeneracy 2 at each 
empty site. Jsij is the bilinear interaction, K the quadrupolar interaction, 
and /i the crystal field. In the last few years the disordered BEG model 
has been studied for random values of the Eij = ±1 ]TB| , |TU[. Recently the 



Degenerate BEG (DBEG) ||20|, ^ has been found suitable to describe the 
martensitic trasformation. 

It may be useful to write the Hamiltonian of Eq. (|I|) in the following way: 

- (3H = '^[J{ei,jSiSj - l)ninj + r]Jninj] + fi^rii, (2) 

where rj = K/J + 1. For 77 = and J = 00, this model has been extensively 
studied in the last few years to study glassy systems and granular materials 
in the disordered case (i.e. when the Eij variables are randomly distributed 
on the lattice ) [6], |^-||29|. This model can be considered as a model of 
particles with an internal degree of freedom {S = ±1) that interact with an 
effective couphng J^EijSiSj — 1) which is zero for spin configurations that 
satisfy the interactions (i.e. EijSiSj = 1) and gives an infinite repulsion, for 
those that do not satisfy the interaction (i.e. EijSiSj 7^ 1). So these last 
configurations are forbidden for J = 00. 

Here we analyze the model for J finite and 1] ^ 0. For rj there is 
an extra interaction between a pair of n.n. particles,while finite values of J 
correspond to softening the hard core potential between the spin variables. 

In particular we present a mean field analysis of the Hamiltonian of Eq. 
(|1|) in the fully frustrated (FF) case on the square lattice. In this case the Eij 
variables are choosen in such a way that every plaquette (i.e. elementary cell 
of the lattice) is frustrated. In other terms every plaquette has an odd number 
of Eij = —1, so that the four spins of the plaquette cannot completely satisfy 
the interactions. In Fig. 1 we show the Villain ||3^ scheme for the 2D FF 



model, highlighting the differences between the A and B sublattices. For this 
FF lattice we have recently ||^ made a mean field analysis of the Frustrated 



Percolation problem ||3^-||43[|. 

In Sec. 3 and Sec. 4 we write down the equations for site magnetizations 
{rriA and m^) and site densities {Da and Db) and these enable us to find the 
critical lines for the order-disorder transitions in our model for the FF case. 
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For K/J > —1 (i.e. r/ > 0) there is a tricritical point which separates 
the critical hne in two branches, respectively characterized by first-order and 
second-order transitions. On the other hand for K/J = —1 (i.e. ?7 = 0) the 
transitions are second-order for any /x. 

Finally we compare the FF behaviour with that of the original 
Ferromagnetic BEG with and without degeneracy. 

2 Mean field analysis 

We will study the model defined by the Hamiltonian (|l|) by evaluating its free 
energy in a mean field approximation. For convenience we will set k = K/ J 

At each site i of the lattice we have to consider the variables Si = ±1 
and rii = 0, 1. For notation purposes it is useful to introduce a new 4-state 
variable Ui such that {ui} = {n^} » {Si} = {1 t, 1 i, T, i} = {1,2,3,4}. 
We can express the old variables in terms of this new variable by means of 
the relations: riiSi = 6u„i - 5z.,,2 and rii = by^^x + by^^i. 

Moreover, using the index r to denote one of the four states of z/j, we can 
define = {5y.^r)^ i-e. the probability that the site i will be found in the 
state Vi = r. Here the angular brackets represent, as usual, the average done 
with the Hamiltonian of Eq. (|I]). 

To obtain the free energy we evaluate first the internal energy of the 
system, which is the expectation value of our Hamiltonian: 

= jJ2e,,{{5.^,x-5.^,2)i5.„i- 5,^,2))+ (3) 

In the MF context we neglect the fluctuations and can simply put 

SO relation (4) implies 

{-pn^J) = J [e^j {p\ - pI) {p{ - pi) + K {p\ + pI) {p{ + pi)] , (5a) 
i-m) = f^{p\+P2)- (5b) 
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The order parameters we will use in the following are the site 
magnetization and the lattice gas particle density Dj expressed by 



rrii = {SiUi) = {S^^^i - (5^.,2) = p{ - P2, (6a) 
A = {ni) = {S.,,i + S,,,2) =p\+pi, (6b) 



from which we have: 



P\ = \(^i + ^i) ^2 = ^ (A - m) ■ (7a) 



These relations and the equivalence condition p\ — p\, together with the 
normalization Yl,t=iPr ~ ^' ii^Pty- 

P3=P4 = ^(1-A). (7b) 



Moreover wc invoke the typical translation invariancc requirement of 
the MF approximation, taken separately on the two sublattices. Then we 
look for a solution in which all the sites of sublattice A (B) have the same 
probabilities, i.e. — \/ i e A and p\.—pf\/i&B. This solution is one 
of the many occurring in the degenerate ground state. 

Using the translation invariance we can write 

-/3Hab = J [ruAmB + k-Da-Db] 

for the expectation value {—j31-Lij) of the partial Hamiltonian relative to any 
AB ferromagnetic bond, i.e. any ferromagnetic bond {ij) such that i & A and 
j E B. A similar relation holds for all the partial Hamiltonians relative to 
any AA ferromagnetic bond. On the other hand, the expectation value of the 

= -1) is 



partial Hamiltonian relative to any BB antiferromagnetic bond 
given by 

-PHbb = J [-ml + kDI] . 

Therefore, for N sites, since the number of A sites and the number of B sites 
are both N/2, the internal energy is 



1 ^ 

-y 



i=l 
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Jz 1 

— [m\ + 2mAmB - m% + k {Da + -0^)^] + -/x {Da + -Db) • 



For the evaluation of the MF entropic term we use the factorization 
property of the probabihty distribution P(z/i, . . . , z/^v) and therefore get 
S^-k E{.} V\nV = -k Eti K In 

Using the translation invariance, this can be written in the form 



S 



kN 2 



-if^iPr^^'Pr+Pr^^Pf)- (9) 



r=l 



Using Eqs. (8) and {^ we can finally write the MF free energy per site 
of the lattice: 

= (10) 



where the probabilities and have to be expressed in terms of the local 
order parameters mA,mB, Da and Db through Eq (7). 



3 Equations for the site Magnetizations and 
Densities 

The knowledge of the free energy allows us to write down easily the MF 
equations that must be satisfied by the order parameters rriA, riiB, Da and 
Db- 

From the stationary relations df/dniA = and df/drriB = it follows 
that 

niA = -DAtanh (^^{^a + ms)^ , uib = -Dstanh (^^{^a - tub) 

' (11) 

Here A = 4J = AJo/kT = T^/T where = 4Jo/^ is the mean field critical 
temperature of the isotropic Ising model recovered by the isotropic version 
of the Hamiltonian (1) in the — > cxd limit. 

Morover from the stationary relations df /ODa = and df /ODb = we 
deduce that 

n2 _ ™2 

kX{Da+Db)+2ii ^ "''A 

a-D.)^' 
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;i2) 



n2 _ ™2 



These relations give in implicit form Da and Db for every and m^. 

Now, replacing Eqs. (|TT]) into Eqs. (|T2|) we get stationarity in the four 
order parameters m^, m^, Da and D^- After straightforward calculations 
we find: 

r, ^ cosh [(A/2)(mA + mg)] 

^ e(-V2){i3A+Ds)-M + cosh [(A/2) (m^ + ms)]' ^ ^ 

sinh[(A/2)(mA + mB)] 

e(-«V2){A4+Ds)-M + cosh [(A/2)(m^ + m^)] ' ^ ^ 

r, ^ cosh [(A/2)(mA-mB)] , . 

B e(-'^V2)pA+Ds)-M + cosh [(A/2) (m^-ms)]' ^ ' 

^ ^ sinh[(A/2)(mA-mi;)] 

B e(-«V2){i3A+Ds)-M + cosh [(A/2) (m^-ms)]' ^ ^ 



These equations can be studied numerically for different values of k in 
order to find the fixed points for every A and ^. This analysis, together with 
the values of the free energy ( [T0| ) for each fixed point, has enabled us to find 
for every /x the critical value Ac where the order parameters , mB , Da 
and Db undergo a first-order or second-order transition. 



4 Critical lines and Results 

We have done our analysis for a number of values of the k parameter, but 
report here, for convenience, only the most interesting cases in the range k ^ 
— 1 (i.e. ?7 ^ 0). Note that the antiquadrupolar phase that generally appears 
in the BEG model for k < does not appear here because our sublattice 
partition is intrinsically different from the usual BEG sublattice partition. 
The critical behaviours are reported in Fig. 2-6 respectively for rj = 1.16, 
1, .84, .5, 0. To appreciate the differences between the FF model and the 
Ferromagnetic model (i.e. eij = 1 for all bonds), each figure contains the 
(a)-section in which we report the behaviour of the Degenerate FF BEG 
model and the (b)-section relative to behaviour of Degenerate Ferromagnetic 
BEG model. In the (a)-section for each T/Tc we give the field — /i/ A were the 
transition from the high- field disordered phase (m^ = ms = and Da = Db) 
to the low-field ordered phase {tjia > 7^ and Da > Db) takes place. 
Bold (dotted) lines represents second-order (first-order) transitions. Dashed 
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lines represent the spinodals, i.e. the boundaries of areas of metastabihty that 
surround any first-order transition hne. Below the first-order transition line, 
the metastable phase is the disordered phase, above this line the metastable 
phase is the ordered phase In the (b)-section for each T/T^ we give the 
field — /i/A were the transition from the high- field disordered phase (m = 
and D < 1/2) to the low-field ordered phase (m > and D > 1/2) takes 
place. As for the (a)-section, bold (dotted) lines represents second-order 
(first-order) transitions; dashed lines represent the spinodals. Fig. 2-6 is 
relative to decreasing values of the extra-interaction t] = k + 1. The overall 
feature is that decreasing rj we obtain a smaller ordered region. This is 
expected if we look at the Hamiltonian (|]) since rj is the extra interaction 
among the particles. In the insert of Fig. 3b and Fig. 4b we report also the 
behaviour of the original BEG. 

For the Ferromagnetic Degenerate BEG we find that the degeneracy 
reduces the area of the ordered region and increases the area of the region of 



first-order transitions, in agreement with recent results [^, ^ 



On the other hand it is known that the frustration has the conflicting 
effect of reducing this region both for the original BEG with random bonds 



T2| and for the DBEG with random field |21|. Here we find that the 
frustration reduces the ordered region and moves the tricritical point toward 
low temperatures, i.e. the frustration in the Fully-Frustrated model (in spite 
of the small degeneracy present) reduces the first order region. 

These results may be useful to study the effects of the softening of 
the hard core potential and the effect of the attraction between particles 
for systems described by Hamiltonian (|^) such as for example glasses and 
granular material. 
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Figure captions 



Figure 1 2d FF model on the square lattice. Straight (wavy) lines represent 
ferromagnetic (antiferromagnetic) interactions, z = 4 ferromagnetic 
interactions start from each site of the sublattice A (open circles); z/2 
ferromagnetic interactions and z/2 antiferromagnetic interactions start 
from each site of the sublattice B (closed circles). 

Figure 2 (a) Critical lines for the FF lattice for k = +.16 (i.e. t] = 1.16). 
Bold (dotted) lines represents second-order (first-order) transitions. 
Dashed lines represent the spinodals. (b) Corresponding critical lines 
for the ferromagnetic model (i.e. Eij = 1 for all bonds). 

Figure 3 (a) Critical lines for the FF lattice for k = {rj = 1). 
The tricritical point is located at T/Tc ~ 0,233 and — /u/A = 
(1/A) In (— 1 + A/\/2) ~ .166. (b) Corresponding critical lines for 
the ferromagnetic model. The insert reports the critical lines for the 
original BEG 0. 

Figure 4 (a) Critical lines for the FF lattice for n = —.16 [rj = .84)., (b) 
Corresponding critical lines for the ferromagnetic model Ferromagnetic. 
In the insert we report the corresponding critical lines for the original 
BEG 0. 

Figure 5 (a) Critical lines for the FF lattice for k = —.5 [rj = +.5), (b) 
Corresponding critical lines for the ferromagnetic model. Observe 
that both in the ferromagnetic and fully-frustrated case the first- 
order transition line continues in the ordered phase, below the tricrical 
point, similarly to the corresponding behaviour of the original BEG 



Figure 6 (a) Critical lines for the FF lattice for k = —1 [rj = 0), (b) 
Corresponding critical lines for the ferromagnetic model. Observe that 
the first-order transition line now disappears, differently from what 
happens in the spin glass case [118], O. 
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